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Abstract: We calculate single partem distributions inside quark and gluon jets within the coherent 
branching formalism, which resums leading and next-to-leading logarithmic contributions. This for- 
malism is at the basis of the modified leading logarithmic approximation (MLLA), and it conserves 
energy exactly. For a wide preasymptotic range of the evolution variable Y = In [E6/Q Q ], we find 
marked differences in the shape and norm of single parton distributions calculated in the MLLA or 
in the coherent branching formalism, respectively. For asymptotically large values Y > 5 — 10, the 
.£h ! difference in norm persists, while differences in shape disappear. In this way, our numerical study 
delineates the jet energy scale needed for a reliable application of both approaches. We also study 
the dependence of the single parton distributions on the hadronization scale Q and on Aq CD , and 
we calculate within the coherent branching formalism the identified quark and gluon distributions 
inside quark and gluon jets. 
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1. Introduction 

Since the early days of QCD, it has been known that destructive interference between soft gluon 
emissions within a jet suppresses hadron production at small values of the momentum fraction 
x = Ph/E. For inclusive single-parton distributions, this is seen already in the double logarithmic 
approximation (DLA), which predicts a hump-backed plateau [1]. However, if one considers the 
kinematic regime of sufficiently small momentum fractions x, and sufficiently large jet energies 
E, where In [1/x] ~ In [E/Q ] ~ O (l/^/cCj, then corrections to the asymptotic DLA result are of 
relative order ^Ja~ s for the peak position of the hump-backed plateau [1-3] . These corrections remain 
sizable up to the highest experimentally accessible jet energies. The coherent parton branching 
formalism [4, 5] leads to evolution equations for the inclusive single- and multi-parton intra-jet 
distributions, which contain the complete set of next-to-leading O{y/ol) corrections, as well as a 
subset of higher order corrections. The modified leading logarithmic approach (MLLA) [5, 6] to 
inclusive parton distributions, which is at the basis of many phenomenological comparisons, can be 
obtained from these evolution equations after further approximations. 

One may ask whether the evolution equations of the coherent branching formalism contain more 
physics than MLLA. This idea is not supported by parametric considerations, since both MLLA 
and the original evolution equations are complete up to the same order in JoT s . However, the idea 
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is supported by kinematic considerations, since the original evolution equations conserve energy 
exactly, while MLLA does not. This has motivated in recent years several works which aim at going 
beyond MLLA. In particular, in an approach referred to as NMLLA [7-9], one keeps all terms of 
the original evolution equations, which are one order 0{^/a~ t ) higher than the MLLA accuracy. In 
the same manner also 0(a s ) corrections to MLLA have been calculated [10,11]. On the other hand, 
one may solve the original evolution equations numerically without any further approximation. So 
far, this has been done for fully integrated partonic jet multiplicities only [12,13]. The main result 
of the present work is to solve these original evolution equations for the inclusive single-parton 
distributions, and to compare the numerical results to those of MLLA. 



In the remainder of this introduction, we recall two issues which come up if one compares data 
to the QCD calculations reviewed above. First, the range of applicability of the evolution equations 
is limited to sufficiently high energies, and a full numerical solution will allow us to be quantitative 
about the scale above which a comparison of the formalism with data may be regarded as being 
reliable. Second, it has been emphasized repeatedly that the QCD prediction for the hump-backed 
plateau of the single inclusive distribution remains unaffected by the non-perturbative hadronization 
process as long as hadronization is sufficiently local so that it does not alter significantly the shapes 
of partonic momentum distributions. In particular, the phenomenologically successful comparison 
of the (rescaled) partonic MLLA prediction with inclusive hadronic distributions has been viewed as 
support for a local parton-hadron duality (LPHD) underlying the hadronization mechanism [14,15]. 
On the other hand, hadronization models employed e.g. in Monte Carlo event generators, such as the 
Lund string fragmentation or cluster hadronization lead to a significant further softening of single- 
inclusive distributions during hadronization [16-18]. Any difference between MLLA, NMLLA and 
the exact numerical solution to the evolution equations in the coherent branching formalism may 
point to the possible relevance of further higher-order and/or non-perturbative effects, and may 
thus add to our picture of the hadronization process. 



The paper is organized as follows. In Sec. |2|, we introduce the equations of the coherent branch- 
ing formalism [4, 5] and we define all quantities of interest. In particular, we discuss in Sec. [2.1| 
the general ansatz for the initial conditions of these equations. The MLLA approach is briefly 
recalled in Sec. 172. In Sec. 2~ 3 we describe how we solve the full evolution equations of the coherent 
branching formalism. The numerical solutions are discussed in Sec. |[ In Sec. [3.1| , we present the 
single-particle spectra in the phenomenologically accessible as well as in the asymptotic range of 
jet energies. In Sec. |372| we show the results for total multiplicity and total energy of partons. 
In Sec. |3l| we study the dependence of these results on the hadronization scale. In Sec. j3^H we 
explore finally to what extent characteristic differences between the coherent branching formalism 
and MLLA remain, if one allows for the freedom of adapting the norm and hadronization scale 
independently in both approaches. Our method of solving the equation of the coherent branching 
formalism enables us also to determine separately the distributions of quarks and gluons in a quark 
or a gluon jet. The corresponding results are discussed in Sec. f|. Our conclusions are summarized 
in Sec. |[ 
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2. Evolution equations for single parton distributions at small x 



We want to calculate the single distributions D q (x, Y) and D g (x, Y) of partons in a quark or gluon 
jet, respectively. Here, x = p/E is the jet momentum fraction of partons, and the variable 



Y = ln Q~ 



(2.1) 



is written in terms of the sufficiently small jet opening angle 9, and the hadronic scale Qq. We 
denote the logarithm of the momentum fraction by 



I = In - . 

x 



(2.2) 



In the coherent parton branching formalism, the evolution equations for these distributions are 
given by [4, 5] 



dyD,(x,Y) 



7T 



Dq (J Y + In z) + D g , Y + ln(l - *)) - D q (x, Y) 

D g g, Y + In z) + D q (-^, F + ln(l - *)) - D s (x, F) | , 



7T 



P> ff (|,y + lnz) + D g [^—,Y + Hl -z)) -D g (x,Y) 



+ 2n f P qg {z) 



D q (J y + In *) + A, , y + ln(l - *)) - D g (x, Y) 



(2.3) 



with the unregularized splitting functions 

1 + z 2 



P g Jz) = 2C A 



1 — z z 



+ 



+ z(l-z) ) , 

1 



z 1 — z 
P qg (z) = T R (z 2 + (1 - z) 2 ) , T R = -, 

^gg^) = CV ^7 — ■ 



(2.4) 
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Figure 1: The range of z- values, which are kinematically allowed in the evolution of the single-parton 
distributions D q (x,Y), D g (x,Y) up to Y = 2 (left), Y = 5 (middle), Y = 10 (right). 



To calculate jet multiplicity distributions with next-to-leading logarithmic accuracy, one can use 
a DGLAP chain of 1 — > 2 parton branchings which follows an exact angular ordering prescription (in 
contrast to the strong angular ordering prescription in DLA) and in which the coupling constant 
depends on k\ at each vertex. In the set of evolution equations (2^3) the angular ordering is 
implemented by the choice of the arguments (Y+hiz) and (Y+hi(l — z)) of the parton distributions 
and also the k± prescription is applied for the coupling [4]. 

The limits Z-, z + on the z- integral in (|2.3|) are set by the requirement that the transverse 



momentum k± is sufficiently large for perturbative evolution to be valid. That means, k± is larger 
than a hadronic scale Qo 



k± « z(l - z)E6 > Q c 



K + lnz + ln(l -z) > 0. 



(2.5) 



For Y < In 4, this inequality is not satisfied for any real z, so evolution occurs only for Y > In 4, 
where the inequality is fulfilled for 



= l -(l-^l-Ae-y^ <z<z+ = l(l + ^/l-4e- Y y 



(2.6) 



This kinematical regime is depicted in Fig. |T|. One can check that if ( |2.6| ) is satisfied, then the 
arguments (Y + hxz) and (Y + ln(l — z) ) of the functions D q and D g in ( |2.3[ ) cannot be negative. 
The coupling constant can be written as 



a s (k 2 j 



7T 



2N C (3 Y + \nz + \n(l- z) + \' 



^-N c - —rifTji 



A = In 



Qo 



A 



QCD 



(2.7) 



2.1 Initial conditions and ansatz for solution 



A parton shower in the real world is pictured often as a quark or a gluon produced initially with 
high energy and virtuality and evolving perturbatively from large Y to a small hadronic scale 
Y . Evolution occurs by emitting partons at smaller and smaller angles, till a minimal angle 8 Q is 
reached, at which non-perturbative hadronic effects set in. The numerical solution of the evolution 
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equations ( |2.3| ) will proceed in the opposite direction, that is, the initial conditions for the solution 
of ( |2.3| ) are set at the small scale Yq which corresponds to the final partonic state of the physical 
process. These initial conditions specify which partons are measured. For instance, consider the 
initial condition 

D q (x,Y ) = 6(1 - x) , D g (x,Y )=0. (Case I) (2.8) 

This initial condition for the evolution equation specifies that the partons at the low scale Yq cannot 
split further and, by virtue of the choice ( |2.8|) , these partons are quarks. Starting the evolution from 
the initial condition (|2.8| ), the functions D q (x, Y) and D g (x, Y) denote the single- quark distributions 
in a quark and a gluon jet, respectively. Alternatively, for the initial condition 

D q (x,Y Q ) = 0, D g (x,Y Q ) = 6(l-x), (Case II) (2.9) 

the functions D q (x,Y) and D g (x,Y) denote the single- gluon distributions in a quark and a gluon 
jet, respectively. Previous studies [5-13] have focused on the distribution of all partons in a quark 
or a gluon jet. This is obtained by evolving D q (x, Y) and D g (x, Y) from the initial condition 

D q (x, Y ) = 8(1 - x) , D g (x,Y ) = 5(1 - x) . (Case III) (2.10) 

In Sec. [5] of this paper we present the results for the case of the initial condition (III) whereas in 
Sec. £| we consider cases (I) and (II). 

The most general ansatz for the solutions of the evolution equations ( |2.3| ) reads 

D q (l,Y) = S q {Y)S(l)+Q(l,Y), (2.11) 
D g (l,Y) = S g (Y)5(l) + G(l,Y). (2.12) 

Here, the discrete parts are expressed in terms of the Sudakov form factors S q (Y) and S g (Y), which 
denote the probability that a parton does not split in the evolution between Yq and Y. These 
discrete terms arise as a natural consequence of the initial conditions (|2.8j) , ( |2.9|) and ( |2.10| ). The 



continuous parts Q(l,Y) and G(l,Y) vanish at the initial scale Y = Yq and arise from the further 
evolution with fl2~3[). 

2.2 The MLLA evolution equations 



In the following sections, we shall compare the numerical solution of (|2.3|) to the results obtained 
in the MLLA approach. Here, we recall the main approximations involved in deriving MLLA from 
the Eqs. Q2.3|) . For more details about MLLA, we refer to [1,5]. 

In the limit of large Y and small x one can perform an expansion of Eqs. ( [2.3|) in powers of 
y/a^. Formally, this is done by treating logarithms Y and In 1/x as being of order 1/ y/oT s . Keeping 
only terms of order O(yfal) on the right hand side of Eqs. (|2.3|) (i.e. keeping terms of order 0(1) in 
the distributions D q and D g ), these equations reduce to the DLA evolution equations with running 
coupling. The subleading corrections to DLA are obtained by keeping terms up to 0(a s ) on the 
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right hand side of (|2.3|) (thus keeping terms of order O(^fal) in D q and D g ). This leads to [19] 
d Y D q (x, Y) = ^ { jf * ^ 7o 2 (F + In z) D g (~ Y + In *) - a qg 7o 2 (F) D 3 (*, Y) } , (2.13) 

9 yJ D g (x,r) = / — 7o(^ + ln^) D g (-,Y + In z) - a gql 2 (Y) D q (x,Y) - a ggl 2 (Y) D g (x,Y) , 

J Z \Z J 

(2.14) 



where ^(Y) = l/(/3(Y + A)) and 



3 2n f TR 1 f U M ^ 4 t 1 ,9 -in 



In deriving Eqs. ( p,13| ) and ( |2.14| ), the semi-hard splittings z ~ 1 from the evolution equation (|2[ 
have been taken into account only partially. For this reason, energy is not conserved exactly in the 
above equations. 

Formally, to arrive at the MLLA equation [1,5], one requires the DLA relation 

D q (l,Y) = ^-D g (l,Y). (2.16) 

By inserting this expression into ( |2.14p , one obtains a closed equation for the parton distribution 
in a gluon jet 

-1 dz 



d Y D g (x,Y)= [ — 1 2 (Y + \nz)D g (-,Y + \nz)-a ll 2 (Y)D g (x,Y), (2.17) 
Jo z \z / 



where 



a l — &gg T Cl gq 



3 3 i \ N r 



(2.18) 



Eq. (|2.17|) is commonly referred to as MLLA evolution equation. From this, the parton distribution 
in a quark jet is usually determined with the help of ( |2.16| ). Compared to the DLA approximation, 
the MLLA equation shows mainly two improvements: the coupling constant is running and the 
negative term oc — ■y 2 l (Y)G (x, Y) on the right hand side of Eq. Q2.17D accounts for recoil effects 
which lead to a softening of the spectra as compared to DLA. 

We note that the MLLA evolution equation fl2.17|) propagates only the initial conditions (II) or 
(III), which give the same result for D g (x, Y). This is because within the MLLA approach one cannot 
differentiate between quarks and gluons in a quark or a gluon jet. In the following, all numerical 
results for the MLLA evolution are obtained with the initial conditions D g (x, Y) = 5(1 — x). Using 
the ansatz (|2.12|) one obtains from Eq. (|2.17| ) the MLLA Sudakov form factor 

Smlla(Y) = ( — ±-J (2-19) 
and the equation for the continuous part G(l, Y). 
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2.3 Evolution equations in the coherent branching formalism 



In this subsection, we return to the full evolution equations ( |2.3| ) of the coherent branching formal- 
ism. Inserting the ansatz (|2.11| ), fl2.12j ) into (|2.3|) , we find for the Sudakov form factors 



d Y S q (Y) 
dyS g {Y) 



- ~ \ Z+ dz^S± [P qq (z) + P gq (z)]S q (Y), 

2 .L it 



~ [ Z+ dz "iffil \P gg {z) + 2n f P qg (z)) S g (Y). 

2 Jz- 71 



(2.20) 
(2.21) 



The above equations are not coupled and can be easily solved. It follows from the integration 
boundaries z_, z + given in Q2.6Q , that 



S q {Y) = S g (Y) = 1 for Y < In 4 . 



(2.22) 



After exploiting the z 1 — z symmetry of the transverse momentum ( |2.5|) and the integration 
boundaries, and with the help of the property P qq (z) = P gq (l — z), we find for the initial condi- 
tion ( P.10|) and for Y > In 4 



S q (Y) 
S 9 (Y) 



exp 



exp 



C F * 



1 



dY' 



dz 



1 + z 2 



Y' + lnz(l - z) + A l-z 

dz 



(2.23) 



Y r z+ 

dY' / dz 

4AT C /Un4 A_ Y' + \nz(l-z)+X 



X 



2CV 



l-z 



+ — z +*(!-*))+ 2n f T R {z 2 + (1 - ^) 2 ) 



(2.24) 



We have performed the double integrals in Eqs. ( p.23 ) and (|2.24|) numerically. Fig. § shows the 
result for three values of A. For Y > In 4, these factors decrease almost exponentially in sharp 
contrast to the MLLA behavior ( |2.19|) . One can show that for large Y the slope changes with Y 
according to 



d\n S, 



1,9 



dY 



- -Cq,g In 



Y + X 



(2.25) 



where c q = Cf/{N c [3) and c g = 1/(3. 

The solutions S q (Y), S g (Y) enter the evolution equations for Q{1, Y) and G(l, Y). After chang- 
ing the integration variable z to I = \n(z/x), the equations for the continuous parts of the spectra 
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Figure 2: The Sudakov form factors (|2.23|) and (p.24|) for three different values of A. 



can be written as 



dyQ(l,Y) 



dyG(l,Y) 



dl T s (l -l,Y)< P qq (l - Q(l, Y-l + l) + P gq (l - I) G(l, Y-l + l) 



+ T s (l, Y)ip qq (l) S q (Y - I) + P gq (l) S g (Y - I) | + cVhi S q (Y) Q(l, Y), 



dl T s (l -l,Y)< P gg (l - I) G(l, Y-l + l) + 2n f P qg (l - I) Q(l, Y-l + l) 

i x 

+ T s (l, Y)^P gg (l) S g (Y -l) + 2n f P qg (l) S q (Y - I) j + d Y \n S g (Y) G(l, Y), 



where we have introduced the shorthand 



r s (i,Y) 



1 



2N C (3 Y -l + \n(l-e- l )+X 



The limits of integration in Eq. ( |2.26| ) are 

2 



I -In 



1 - Vl - 4e" 



l-ln- 



1 + Vl - 4e" y 



(2.26) 



(2.27) 



(2.28) 



We recall that the functions G(l,Y) and Q(l,Y) vanish when the first argument is negative. We 
also note that the form (|2.26|) for the evolution equations in the coherent branching formalism is 
valid for all initial conditions and the logarithmic derivatives dyhiS qtg (l,Y) on the right hand side 
of these evolution equations are formal notational shortcuts, denoting the expressions fl2.23|) , ( j2.24| ) 
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irrespective of the initial conditions. These logarithmic derivatives decrease only logarithmically 
with Y, as seen from Eq. ( |2.25| ). Hence, the rapid decrease of the Sudakov form factors itself does 
not render the corresponding terms in ( |2.26| ) negligible at large Y. 



3. Numerical results 

We have solved numerically the evolution equations for the continuous parts of single particle 
distributions in quark and gluon jets in the coherent parton branching formalism ( 2.26| ) and in the 
MLLA approach (|2.17|) . The solutions depend on 

A = lnQ /A Q CD, (3.1) 

which specifies in units of Aqcd the value of the hadronization scale Qo up to which the parton 
shower is evolved. It is a peculiar feature of the MLLA approach, that its solutions have a well- 
defined finite limit A — > 0, which is known as the "limiting spectrum" and which is at the basis of 
many phenomenological comparisons [1,14,15,20,21]. In general, however, one cannot expect to 
obtain finite results in the limit A — > 0, since A regulates the Landau pole of the strong coupling 
constant ( |2.7|) . In fact, we shall find that the solutions for the evolution equations ( |2.3|) can be 
obtained for a finite value of A, only. 

In the following, we show numerical results for A = 0.01, A = 0.2 and A = 1.0. The MLLA 
limiting spectrum is approximately 3 — 4% larger in norm and similar in shape, compared to the 
MLLA solution for A = 0.01. Hence, our choice for this smallest value of A is motivated by the 
interest in exploring the solutions of the coherent branching formalism (|2.3|) for a parameter range 



in which the MLLA limiting spectrum is almost reached. Our choice of the largest value A = 1.0 is 
motivated by the fact that parton showers implemented in Monte Carlo event generators typically 
end the perturbative evolution at a hadronization scale Qo which is significantly larger than Aqcd 
and for which A ~ 1 is an order of magnitude estimate. We have chosen one value A = 0.2 in 
between, and it is for this value that we shall first present our results. 

3.1 The distributions Q(l,Y) and G(l,Y) 

In Figs. |H and |], we present results for the continuous part of the single-parton distributions inside 
a quark and a gluon jet, evolved from the initial condition ( |2.10| ) at Yo = up to different values of 
Y. Fig. ^ shows results for the evolution up to values of Y < 5 whereas Fig. f| focuses on evolution 
to much larger values 15 < Y < 30. The relation between the evolution variable Y and the jet 



energy E depends on the jet opening angle 9 and the hadronization scale Qo, as seen from Eq. ( |2.1| ). 
In the case of large angles, the approximate formula ( |2.1|) must be replaced by the exact relation 
Y = hx(Esind/Qo)- For a given jet energy, the maximal value of Y is obtained with the maximal 
choice for 9 and the minimal choice for Qo- For instance, for a maximal angle, 9 = n/2, and a very 
low value Q = 200 MeV, the value of Y = 5 corresponds to an energy of E = 30 GeV and Y = 7 
corresponds to E — 220 GeV. On the other hand, for a smaller opening 9 = 0.3, a choice of Y = 5 
corresponds to E — 100 GeV. Also, for 9 = 0.3, Y = 5 corresponds to E — 250 GeV if one chooses 
a larger hadronization scale Qo = 500 MeV. These numerical examples illustrate that while the 
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Figure 3: The continuous parts of the inclusive distribution of all partons in gluon (upper panel) and 
quark jet (lower panel). Results are shown for evolution in the coherent parton branching formalism (p.3[), 



and in the MLLA approach (2.17) for A = 0.2 and Y = 2.5, 3.3, 4.2 and 5.0. 



correspondence between jet energy and Y may vary significantly depending on the choice of Qq, the 
range of Y < 5 — 7 matches quite well the entire range of jet energies which have been compared so 
far to data on inclusive single hadronic distributions inside jets. This is one reason for considering 
first the evolution up to Y = 5. We note that similar numerical estimates relate Y = 10 to a minimal 
jet energy E > 4 TeV. So, our choice of larger values of Y can be motivated solely by the theoretical 
interest of testing the asymptotic behavior of these evolution equations. In addition, we recall that 
the evolution equations (|2.3|) have been derived in an approximation, in which logarithms are large 
and corrections of order Y ~ \j ^foT s are kept. Starting at Yq, some evolution will be necessary 
to satisfy this condition. As we shall see, this initial stage of the evolution leads to some peculiar 
features in the solution which can persist up to values of Y = 5 or larger. 

We consider first the case for Y < 5. As seen in Fig. |^, the single inclusive distributions 
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I = In 1/x 

Figure 4: Same as Fig. [|, but for evolution up to Y = 15, 20, 25 and 30. 



obtained for the evolution within the coherent branching formalism show for small Y a large yield 
at low values of I <C Y. This effect is more pronounced for the distribution in a quark jet, than 
for the distribution in a gluon jet. It can be understood on general grounds from the fact that 
the initial conditions ( |2.10| ) for the evolution are delta-functions at I = 0, and that the results for 
small Y remember these initial conditions. Moreover, quarks are less likely to branch than gluons, 
and in their first branching they are likely to keep a large fraction of their energy. This tends to 
make the distribution in a quark jet harder than in a gluon jet, as clearly seen in Fig. |||. Upon 
further evolution to Y — 3.3, 4.2, 5.0, the yield of hard partons (say partons with I < 1) in a gluon 
jet decreases rapidly to a value comparable with (and for some region of I even smaller than) the 
results obtained in the MLLA approach. On the other hand, the yield of hard partons in a quark 
jet remains enhanced significantly over the yield obtained in the MLLA approach up to relatively 
large values of the evolution parameter Y. We recall that the parametric reasoning underlying the 
evolution equations (|2.3| ) is valid for sufficiently large Y and I ~ 0(Y). However, also in the region 



- 11 - 




5 10 15 20 25 30 

Y 



Figure 5: The evolution in Y of the peak position of the continuous parts of the inclusive parton distri- 
butions. Results are shown for all partons in a gluon jet (upper straight line) or a quark jet (lower straight 
line), calculated in the coherent parton branching formalism ( |2.3| ). Results are also shown for the MLLA 
approach fl2.17| ) (dashed line), for which the peak positions of the distributions for quark and gluon jets 
coincide. All results are calculated for A = 0.2 and the peak position Z max is plotted in units of Y. 



around the peak of the single inclusive distributions (which will lie for sufficiently large Y within 
this kinematic range of validity), there are significant differences between the MLLA approach and 
the coherent branching formalism. 

From Fig. we have seen several preasymptotic features, which become less pronounced upon 
further evolution in Y. Fig. |] illustrates the degree to which these preasymptotic features disappear 
with increasing Y. The numerical results for the coherent branching formalism and the MLLA ap- 
proach agree very well in shape and differ in norm by a factor of order unity which is ^-independent 
to good accuracy. Closer investigation shows that the factors needed to adjust the norm between 
both approaches differ for gluon jets (normalization factor 1.24 for A = 0.2) and quark jets (nor- 
malization factor 1.30 for A = 0.2). Fig. |4] illustrates the statement that in the asymptotic region of 
evolution towards sufficiently large Y, the region around the peak of the single-parton distributions 
lies at sufficiently large I to be calculable within the accuracy of the single branching formalism. 
In this asymptotic region, MLLA and coherent branching formalism agree up to the norm. On 
the other hand, Fig. |3] illustrates that for A = 0.2, this asymptotic region is not yet reached for 
Y < 5 — 7, which is the range in which the bulk of the currently known experimental data lie. 

To further illustrate the slow but steady convergence in shape between the MLLA results for 
single-parton distributions and those of the coherent branching formalism, we have plotted in Fig. |5] 
the evolution of the peak position of these distributions as a function of Y. We see that the peak 
positions in both approaches differ largely for Y < 5 — 7, consistent with our earlier remarks. 
In general, the gluon jets split more than quark jets and thus show a softer distribution and a 
peak at larger value of / max . This difference, however, decreases with increasing Y, as does the 
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difference with the peak position in the MLLA approach. The value of l mauX /Y decreases slowly 
with Y, but even for the largest values explored here (i.e. Y = 30 corresponding to the jet energy 
of E ~ 10 12 GeV), the position of the peak differs significantly from the DLA result l max /Y = 1/2. 

3.2 Total multiplicity and energy conservation 

The single-parton distributions inside quark and gluon jets can be characterized by their moments. 
The zeroth and first moments of the inclusive single-parton distributions are of particular interest, 
since they define the total parton multiplicity in a quark or gluon jet, 



N M (Y) = JdlD q>g (l,Y), (3.2) 
and the total energy fraction inside a quark or gluon jet 

i E p( Y ) = JdxxD qj9 (x,Y). (3.3) 

partons 



For an energy- conserving evolution of D q ^ g , the first moment of the single-parton distribution must 
yield the total jet energy E, and the fraction (|3.3p must equal one. The bottom panel of Fig. 



illustrates, that this is satisfied for the coherent branching formalism. Since we know that the 
evolution equations ( |2.3p conserve energy, this is a cross check of the numerical accuracy of our 
solution. We note that the value of the first moment ( |3.3| ) is determined largely by D q ^ g (x, Y) in the 
region of large x. However, parametrically, the evolution equations ( |2.3|) have been derived for the 
region of small x, where the logarithm I is large. Fig. ^ also illustrates that the MLLA approach 
does not conserve energy exactly. So, on kinematic grounds, the coherent branching formalism is 
certainly preferable. We note, however, that for large values of Y, as seen in Fig. |], the region 
around the peak of the distribution E w will be at x values, which are much softer than the region 
of x space which contributes dominantly to ( |3.3| ). For this reason, a formalism which does not 
conserve energy may account accurately for the functions D q g (x, Y) in a wide region around the 
peak, if Y is sufficiently large. 

Within the MLLA approach, one customarily relates by a prefactor the single-parton distribu- 
tions in quark and gluon jets, so that G(x, Y) = N c /Cp Q(x,Y). This relation is valid to DLA 
accuracy and leads to the first moment ( |3.3j ) of a quark jet equal to the total energy faction of a 
gluon jet multiplied by 4/9. We did not include this result in Fig. |6|, since it does not add further 
information. 

We now turn to the total parton multiplicities within a quark or a gluon jet, plotted in the upper 
panel of Fig. These multiplicities satisfy evolution equations for integrated parton distributions, 
which are much simpler than the set of equations ( |2.3| ), and which have been studied numerically [12, 
13]. We checked that we reproduce these results, if we adopt the same modifications to the coherent 
branching formalism (|2.3j ). We also tested the numerical accuracy of our routines for D q ^ g (x, Y), by 
comparing the integral ( |3.2|) to the result of the simpler evolution equations for multiplicities, and 
we found perfect agreement. Within MLLA, the multiplicities for the limiting case (A = 0) and at 
large values of Y are known to rise like [1,2,7] 

Nmlla(Y) ~ Y - a ^ +l ' 4 e 2 v 7 ^ . (3.4) 
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Figure 6: Upper plot: Independence of the total parton multiplicity ( |3.2| ) in a quark or gluon jet, 
calculated in the coherent parton branching formalism ( |2.3D (straight lines), and in the MLLA approach 
( 2.17] ) (dashed lines) for A = 0.2. Lower plot: The ratio ( |3,3| ) of the total jet energy contained in the single 
parton distribution to the initial jet energy. The upper panel is for gluon jets, the lower one for quark jets, 
A = 0.2. 



This strong rise of the multiplicity is seen in Fig. ^. At large Y, the rise is the same in the MLLA 
approach and in the coherent branching formalism. In comparison to the MLLA approach, the 
multiplicities obtained in the coherent branching formalism are a factor ~ 1.24 for gluon jets and a 
factor ~ 1.30 higher for quark jets at large Y. This is consistent with our observation in Fig. f| and 
indicates that at large Y, the distributions obtained in both approaches can be made to coincide 
by an approximately ^-independent multiplicative factor. 
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Figure 7: Same as Fig. but for A = 0.01. 



3.3 Dependence of Q(x, Y) and G(x, Y) on A 



The numerical results shown in the previous subsections |3.1| and |3.2| were obtained for A = 0.2. 
Here, we discuss the dependence of these results on A. From the point of view of data comparison, 
it is natural to regard the jet energy E and the hadronization scale Q as the two independent 
variables, while Aqcd is fixed. Specifying the F-range of the evolution and the value of A amounts 
to specifying the variables E and Qq. On the other hand, regarding the values for the jet energy and 
the F-range of the evolution as given, decreasing A amounts to increasing Aqcd- As a consequence, 
a smaller value of A is expected to lead to a more 'violent' evolution (i.e. more parton branchings) 
and it will thus lead to a higher jet multiplicity within the same y-interval. 

This general expectation is confirmed by the curves in Figs. [F] and [8], which show the distribu- 
tions of all partons in a quark or a gluon jet for evolution over small ^-intervals up to Y < 5. In 
particular, one sees that after evolving for a few units in Y, the parton distributions for A = 1.0 are 
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Figure 8: Same as Fig. but for A = 1.0. 



still peaked at very large values of x, where In 1/x <C 1. Qualitatively, the structures at large x are 
due to the fact, that for small Y the phase space for evolution shown in Fig. [I] excludes evolution 
to very small values x, and the driving term of the evolution is the Sudakov form factor at x — 1. 
Upon evolution in Y, this Sudakov form factor decreases faster for the more violent evolution with 
A = 0.01 than for less violent one with A = 1.0, as can be seen from Fig. |2|. For this reason, the 
parton distributions in Fig. |7|, which were calculated with A = 0.01, decrease under evolution in Y 
much faster in the region of larger x, than distributions evolved with larger values of A, shown in 
Figs. [I] and ||. The same qualitative features (though quantitatively less pronounced) are found in 
the MLLA approach, where for small Y the continuous part of the single parton distribution differs 
significantly from zero at large x, in particular for the case A = 1.0 of a less violent evolution. 

For A = 0.01 and A = 1.0, we have calculated the single-parton distributions inside quark and 
gluon jets for evolution up to Y = 30 (results not shown). In this asymptotic regime, similarly 
to the case A = 0.2, the results obtained from the coherent branching formalism and the MLLA 
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approach can be made to coincide by a simple adjustment of the norm. For A = 1.0, the MLLA 
results are at large values of Y a factor 1.11 smaller for gluon jets and a factor 1.16 smaller for 
quark jets. For A = 0.01 the corresponding factors are 1.72 for gluon jets and 1.80 for quark jets. 
As noted before in studies of total jet multiplicities [12,13], the coherent branching formalism does 
not allow for a limiting curve, A — > 0, while the MLLA approach does. Our numerical findings are 
consistent with this observation and indicate that for decreasing A, the difference in norm between 
both approaches increases slowly but without bound. 

For a more violent evolution, i.e. for smaller A, single parton distributions increase in multi- 
plicity and become softer. As a consequence, the peak position of the distributions is expected to 
shift to larger values of In 1/x as A decreases. This expectation is confirmed in Fig. |9], which shows 
the peak positions of single-parton distributions for A = 0.01 and A = 1.0. 



3.4 Matching MLLA to the coherent branching formalism 

As shown in Figs. |3], [7| and |8|, single parton distributions calculated in the MLLA approach differ 
for Y < 5 both in norm and in shape from those calculated in the coherent branching formalism. 
Besides the normalization, the spectra depend on two additional parameters, namely Qo and Aqcd 
or equivalently Y and A. Here we ask to what extent a variation of the three parameters is sufficient 
to make results from the coherent branching formalism coincide with a given result of the limiting 
fragmentation function (A — > 0) obtained in the MLLA approach. One motivation for this study 
comes from the observation that the MLLA limiting fragmentation function is at the basis of many 
phenomenologically successful comparisons. So, even if the comparison of QCD predictions for 
single-hadron distributions in jets is not restricted to the calculation of single-parton distributions 
and involves difficult issues in the modeling of hadronization, one wonders to what extent the shape 
of partonic distributions in the two calculational schemes can be made to agree by a suitable choice 
of parameters. 



Fig. [XO shows an example of the extent to which this is possible. Adjusting the norm, the 



evolution time, Y, and the value of A, the single-parton distributions in a quark jet can be made 
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Figure 11: The ratio G(l,Y)/Q(l,Y) in the coherent branching formalism for two different values of A 
and evolution up to different Y = 15, 20, 25 and 30, from top to bottom. By construction, MLLA results 
assume the DLA relation fl2.16| ). 



to coincide almost for both formalisms at Y = 5. More generally, one sees from Figs. [| and |9] that 
the peak position of the single-parton distributions in either a quark or a gluon jet can always be 
adjusted to agree between both approaches by varying Y and A in the coherent branching formalism. 
The total yield can then be adjusted by the norm. 

On the other hand, it is not possible for Y ~ 5 and with the same choice of A to get results of 
both formalisms to coincide for the distributions in both quark and gluon jets . This is illustrated 
by the example given in Fig. [HJ where marked differences persist between the distributions in 
gluon jets, once the distributions in quark jets have been adjusted. This conclusion is supported 
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more generally by Fig. |TT], which shows the ratio of single-parton distributions of a gluon and a 
quark jet. Within the MLLA approach, these distributions vary in norm but not in shape, since 
the ratio G(l,Y)/Q(l,Y) is fixed to the DLA value N c /Cf- On the contrary, within the coherent 
branching formalism, the distributions in a quark jet and a gluon jet vary both in norm and shape. 
In this latter case, the deviation from N c /Cp is more pronounced in the region of relatively large 
x, and upon evolution in Y this region becomes disentangled from the region around the peak of 
the distribution. In short, while there are marked differences between both formalisms at all Y, for 
sufficiently large Y these differences die out in the region of large I ~ 0(Y) to which the coherent 
branching formalism and the MLLA approach are tailored. 

4. Identified parton distributions in quark and gluon jet 

So far, we have discussed the distribution of a//partons within a quark or a gluon jet. In the coherent 
branching formalism, we can study separately the distribution of quarks and gluons within a quark 
or a gluon jet. We characterize these distributions by subscripts. For instance, Q g (l,Y) is the 



single-gluon distribution in a quark jet. As explained in subsection |2.1| , single-quark distributions 
are obtained by evolving from the initial conditions ( |2.8| ) and single-gluon distributions are obtained 
by evolving from the initial condition Q2.9|) . The phase space constraints imply that there is no 



evolution in the coherent branching formalism up to Y = In 4. Since the continuous parts of the 



distributions vanish by definition at Y — In 4, the initial conditions (|2.8| ) and ( |2.9|) translate directly 
to the initial conditions for the Sudakov form factors 

S g (Y — In 4) = , Sq(Y — In 4) = 1 , if only quarks are counted, ,^ . 

S g (Y — In 4) = 1 , S q (Y — In 4) = , if only gluons are counted. 

It is straightforward to check that if a Sudakov factor equals zero at In 4, then it stays zero at any 
value of Y. 

In general, the distribution of quarks inside a jet is harder than that of gluons. This is a 



consequence of the fact that gluons are more likely to split. In the left panel of Fig. [12], which 
shows the identified parton spectra in jets for A = 0.2, this is clearly seen for all values of Y. 
A particularly pronounced feature is seen in the single-quark distribution for quark jets: up to 
Y = 5, this distribution peaks at very large x values corresponding to I < 0.5. The integrated 
total quark multiplicity in this distribution, shown in the right panel of Fig. [T2|, grows slowly but 
remains of order unity during this evolution, indicating that this single quark distribution follows 
closely the momentum distribution of the evolved parent quark. This parent quark evolves via 
q — > qg splittings predominantly by emitting soft gluons. For this reason, the quark distribution 
remains hard and the total multiplicity in the quark jet is largely dominated by subsequent g —>■ g g 
splittings. It is in precisely this sense that the enhancement at large-x is a remnant of the 5-function 
in the initial condition, which gradually becomes negligible in the evolution to larger Y. 

We have emphasized repeatedly that the approximations involved in deriving the coherent 
branching formalism do not guarantee an accurate description of the large-x region. We note, 
however, that the reasons which we have identified for the pronounced large-x enhancement of 



Q q (x, Y) in Fig. 12 do not refer to particular features of an 0(y/al) approximation, but emerge 
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Figure 12: Left hand side: The continuous part of the single-parton distributions inside a gluon jet 
(upper panel) and a quark jet (lower panel). Different curves denote the single-quark (dashed) and single- 
gluon (dotted) distributions inside a jet, the explicit sum of these two distributions (dash-dotted), and 
its comparison to the distribution of all partons in a jet (solid) obtained by evolving the initial condition 
(fTIC| ). The two latter curves coincide, as expected. Right hand side: the Independence of the total 
multiplicity of quarks, gluons and all partons calculated for the cases shown on the left hand side. 



solely from generic features of the branching of quarks and gluons. As a consequence, we expect that 
these features, though observed outside the strict region of validity of the formalism explored here, 
correspond to physical reality and persist in a more complete formulation of the problem, which 
is accurate in the large- X region. Indeed, in Monte Carlo event generators, such an enhancement 
is seen on the level of partonic distributions [16-18]. In these event generators, it is generally the 
hadronization mechanism which connects the color of the leading quark to the rest of the event, and 
which thus leads to a hadronic single-inclusive distribution which is much softer than the partonic 
one. 

For the distribution of quarks within a gluon jet, the situation is clearly different since the 
splitting function for g — > qq does not give rise to a logarithmic enhancement. For this reason, 
the total quark multiplicity inside a quark or a gluon jet rises only slowly. For gluon jets as well 
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as for quark jets, the jet multiplicity is ultimately dominated by gluons. The distribution of the 
hardest gluon in the first g — > g g splitting process is likely to contribute significantly to the large 
values which G(l,Y) shows in the large x region (say / < 1.0) for small values of Y. Upon further 
evolution in Y, this structure disappears faster in G(l,Y) than in Q(l,Y), since it is dominated in 
the first case by gluons, which split more readily than quarks. 



5. Conclusion 

We have calculated single parton distributions inside quark and gluon jets by solving the full evolu- 
tion equations in the coherent branching formalism. For sufficiently large jet energies, corresponding 
to Y > 5 — 10, results obtained from the MLLA and from the coherent branching formalism agree 
very well in shape. They differ in this asymptotic regime solely by a factor of order unity in norm, 
and this factor depends weakly on the choice of hadronization scale. For smaller values of the jet 
energy (Y < 5 — 7), however, both formalisms show marked differences in norm and shape. We 
have characterized these differences in detail in sections [3] and 

We recall that, parametrically, both MLLA and the coherent branching formalism are valid for 
sufficiently large Y in a kinematic region in which I ~ 0(Y). As a consequence, the differences 
observed between both approaches for Y < 5 — 7 can be regarded as preasymptotic effects, which 
persist solely outside the strict region of validity of both approaches. If one adopts this view, 
then our study has contributed to delineating the region of quantitatively reliable applicability of 
both approaches, namely Y > 5 — 7. In this context, we note that for Y < 5 — 7, single-parton 
distributions in both approaches cannot be made to coincide by a suitable choice of parameters for 



both quark and gluon jets, while they can be made to coincide for one of them, see section |3~4 . 

On the other hand, the MLLA limiting spectrum has been compared successfully to a large data 
set of single inclusive hadronic distributions for jet energies, which lie in the range Y < 5 — 7. These 
data comparisons required the assumption of local parton hadron duality and, by their success, gave 
support to this assumption. In the present work, we have found that for Y < 5 — 7, distributions 
calculated in the coherent branching formalism show a significant enhancement at large x, which 
is much less pronounced and vanishes much faster with Y in the MLLA approach. Also, this 
enhancement has not been observed in the measured hadronic distributions. 

Although these findings were made in the region of small Y and large x, which lies clearly 
outside the region of quantitative applicability of the formalisms in question, we have argued in 
section [| that the structures found at large x in the coherent branching formalism emerge solely 
from generic properties of the branching of quarks and gluons and should thus persist in a more 
complete formulation of the problem, which is reliable at small Y and large x. We note that the 
marked difference between the shape of partonic distributions calculated in the coherent branching 
formalism and the shape of the measured hadronic distributions does not necessarily invalidate the 
coherent branching formalism in the range of Y ~ 3 — 7. Rather, this difference may be indicative 
of a non-trivial dynamics of hadronization, which is characteristically different from a simple one- 
to-one mapping of partons into hadrons at the scale Q ~ C(Aq CD ). Indeed, since the hardest 
partonic components in a jet are connected to the softer ones by color flow, and since hadronization 
is ultimately sensitive to the color structure of the event, the dynamics of hadronization may lead 
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to a significant softening of the single inclusive partonic large-re distributions. Such features are 
realized at least in some models of hadronization, such as the Lund string model. If one adopts this 
view, then our study may be regarded as a contribution to the question to what extent the measured 
inclusive hadronic distributions inside jets lend support to a specific picture of hadronization. 
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